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TOPOLOGICAL INTERACTIONS IN A BOLTZMANN-TYPE 

FRAMEWORK 

ADRIEN BLANCHET * 1 & PIERRE DEGOND 2 


Abstract. We consider a finite number of particles characterised by their posi¬ 
tions and velocities. At random times a randomly chosen particle, the follower, 
adopts the velocity of another particle, the leader. The follower chooses its leader 
according to the proximity rank of the latter with respect to the former. We study 
the limit of a system size going to infinity and, under the assumption of propagation 
of chaos, show that the limit equation is akin to the Boltzmann equation. How¬ 
ever, it exhibits a spatial non-locality instead of the classical non-locality in velocity 
space. This result relies on the approximation properties of Bernstein polynomials. 


1. Introduction 

In this paper, we explore collective dynamics driven by rank-based interactions, i.e. 
that’s to say interactions determined by the rank of the agents with respect to certain 
criterion. There are many examples where such interactions take place. In economics 
for instance, it was extensively analysed in [2] that agents are more sensitive to their 
rank compared to others (salary or wealth for example) than their own independant 
cardinal level. To go further, na studies, in an organisation, compensation schemes 
which pay according to an individual’s ordinal rank rather than their output level. 
Such payoff based on the rank approach also appears very naturally in a variety of 
economics applications such as bids, the labour market, portfolio management, the 
oil market, academic production, reputation, etc. 

Evolutionary game theory studies the evolution of strategies/genes transmitted 
through natural selection. In chimpanzees as in cockroaches a group is formed of 
a dominant male, females and lower order males. Only the dominant male is sup¬ 
posed to mate with the females. However, when the dominant male is absent, the 
females also reproduce with other males giving the preference to males in descending 
order [12]. It is also known that the rank of an offspring strongly depends on the 
rank of its mother |T6], so that in the replicator dynamic process the rank increases 
the chance of reproduction. The study of such models requires taking into account 
interactions depending on the rank of the agents. 

In this article we focus on the dynamics of bird flocks. There is a widespread 
literature of flocking models where the birds react to their neighbours as a function of 
the neighbours’ distance from them within the flock. These are the co-called “metric” 
interactions. In this context, dynamics based on alignment [20] . consensus m or 
attraction-repulsion see have been widely studied. However, there has been 
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recent compelling evidence [Tj that interactions within bird flocks are mostly metric 
free, as the birds react primarily with a limited number of their nearest neighbours 
irrespective of the distances between them. This observation has motivated the 
concept of “topological interaction”, which has been widely echoed in the scientific 
literature ISU7UinUl3Ul9]. 

Our goal is to investigate the large size limit of a system of agents interacting 
through topological interactions. Specifically, we consider a leader-follower model [8J, 
i where at random times a randomly chosen bird, the follower, decides to adopt 
the velocity of another bird, its leader, in the flock. The follower chooses its leader 
according to a probability only depending on the proximity rank of the latter with 
respect to the former. If we assume that the probability has a strong cutoff as soon 
as the proximity rank exceeds a certain value, of the order of seven in actual flocks, 
the considered model is akin to the topological interactions of [1]. 

To our knowledge, [6] is the first mathematical work where interaction rules be¬ 
tween agents depending on their rank are considered. The closest to our work is [T5] 
where kinetic and hydrodynamic models for topological interactions have been pro¬ 
posed. However, the considered dynamics is different from ours. In [15], it is supposed 
that an agent’s velocity relaxes towards an average velocity of its neighbours where 
the relative weights of the neighbours depend on their proximity rank to the con¬ 
sidered agent. Therefore, it is a model of Cucker-Smale type Em combined with a 
topological interaction rule. In [15], a mean-field type kinetic model is rigorously 
derived under some regularisation in the large system size limit and a hydrodynamic 
model under a monokinetic closure assumption is proposed. 

Here, the interaction rule is different and, in the large system size limit, leads to 
a Boltzmann type model with an integral operator describing the balance between 
gains and losses due to the interactions rather than a mean-field model where the 
interactions are described through a force field. From the mathematical viewpoint, 
this makes a considerable difference, as an empirical measure approach is not possible. 
Instead, one has to rely on the propagation of chaos property for the solution of the 
master equation. In the present work, propagation of chaos is assumed and its proof 
is defered to future work. Still, under this assumption, the derivation of the kinetic 
equation is not obvious and as we will see, relies on fine approximation properties of 
Bernstein polynomials. 

Indeed, we will realise that the derivation of a kinetic model requires the estimation 
of the probability that given two particles say numbered 1 and 2, the rank of 2 with 
respect to 1 be equal to a given integer j. Then, the interaction probability of 1 
with 2 in this configuration is a function K(j/N), where N is the total number of 
particles and the function K is characteristic of the considered interaction. Thanks 
to an easy combinatorial estimation, the total probability of 1 interacting with 2 
is found as the Bernstein polynomial approximation of K when N is large. Due 
to some cancellations, the first order correction in powers of 1/N of the Bernstein 
polynomial approximation of K is also needed. This correction can be found in the 
literature [18 ] . 

The paper is organised as follows. In Section [2j we present the N particle dynam¬ 
ics and state the main result. In Section El we derive the master equation of the 
process and the equation for the first marginal under the assumption of propagation 
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of chaos. In Section [4] we precisely state our main result, namely that, in the limit 
N — y oo and the assumption of propagation of chaos, the equation for the first mar¬ 
ginal reduces to a kinetic equation of Boltzmann type with spatial nonlinearity. To 
prove this theorem, we use results on Bernstein’s polynomial approximation from the 
literature [T8]. Section [5] offers some considerations on the limit kinetic equation and 
illustrates our discussion with numerical simulations. Finally, a conclusion is drawn 
in Section [6] 


2. The TV-particle dynamics 


Consider a set of N particles. The particle i is characterised by its position Xi G M n 
and its velocity u* G M n where n > 1 is both the spatial and velocity dimension. For 
a given particle i we can order the other particles relatively to their distance to i. 
More precisely, we have the following: 


Definition 1 (Rank). Consider N particles located at x i,..., apv- Consider the i-th 
particle and order the list (|Xj — a;i|) N by increasing order and denote by 
R N (i, j) G {1,..., N — 1} the position of the j-th item in this list. If two indices 
j and j' are such that \xj — xf — \xj> — Xi\, then we choose arbitrarily an ordering 
between these two numbers. We define R N (i,i) = 0. Now, we define the rank of j 
with respect to i as: 


N 


(.hj) = 


R N (fj) 

N- 1 


N -1 

u 

k =1 


N — 1 


We introduce a function K: r G [0,1] ha- K(r) G [0, cx)) such that 


We define 


K N (r) = 


K{r ) dr = 1 
K(r ) 


2^k =1 l N -1 


in order to have for any i G {1,..., N}: 


N 


N -1 




3 = 1 


k =1 


N — 1 


= 1 . 


In this way, for any i G {1,..., N}, the collection ^ where 

t $ = K N (r N (i,j )), 


defines a discrete probability measure on the set {j G {1,..., iV}, j j- i}. 

Consider N particles {(aq(f), vfit )),..., (x^ft), v^ft))} which are subject to the 
following dynamics (previously referred to as the “Choose the Leader” dynamics 

my- 

- The dynamics is a succession of free-flights and collisions. 
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- During free-flight, particles follow straight trajectories 

r Xi = Vi, 

\ Vi = o. 

- At Poisson random times with a rate equals to N, particles undergo the fol¬ 
lowing collisions process: Pick a particle i in {1,..., N} with uniform prob¬ 
ability 1/N and perform a collision, i.e. pick a collision partner j in the set 
{j G {1,..., N}, j 7 ^ i } with probability nfj and perform: 

f ( Xi,Xj ) remains unchanged, 

\ ( Vi,Vj ) is changed into (vj,Vj). 

Since the rank of j with respect to i is an intrinsic property of the positions of 
the pair of particles and does not depend on how they are numbered, we have the 
following properties of the rank: 

Remark 2. Let (aq,..., x^) be a set on N particles. 

(i) The rank r N (i,j), and, hence , is a function of (aq,..., xjy), i.e. 

r N {i,j) =r N (i,j)(x l ,...,x N ) . 

More precisely, we consider the rank r N {i, j) as a function of L°°(M. nN ). 

(ii) The rank is permutation invariant, i.e. for any permutation a G &n where 

denotes the set of permutations of {1,..., N}, we have 

r N (a(T), a(j))(x a{l) ,.. .,x a[N) ) = r N (i,j)(x 1 ,.. .,x N ). 


The aim of this article is to study the limit of this dynamics when the number of 
particles goes to oo. To do so we will assume that the propagation of chaos property 
holds true i.e. 

N 

f iN) (Z ir -- ,Z N ,t) = HfS\z e ,t), VZGI", VtG[0,oo). 

l= l 

Assuming that fff —> f and := J f^ dv —> p — J f dv, then in the limit 
A^ —> oo, we will prove that / is a solution of the kinetic equation: 

df r 

— (x,v) + v- V x f(x,v) = p(x) / f(x',v)K(M p (x,\x -x\)) dx - f(x, v), 
where M p is the partial mass of p and is defined by 

M p (x, s) = / p(x') dx' , 

J x'£B(x,s) 

and where Bfx,s) = {y G M n | \y — x\ < s} is the ball centered at x and of radius 
s > 0. 


Remark 3. The conservation of mass property holds true by Lemma If applied to 


H = K. 
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In the following section, we derive the master equation for this process, Section [3TT1 
and the first marginal equation for indistinguishable particles, Section 13.21 Then, in 


Section 13.31 we derive the master equation under the assumption of propagation of 
chaos. 


3. Master equation and propagation of chaos 


3.1. Master equation. To simplify the notation, when no confusion is possible, we 
will denote x := (aq,..., xn), v := (iq,. .., Vn), Zi := (ay, iy), Z := (Z i,.. . , Zjv) and 
dZ := doq drq ... d xn d vn- 

As the collisions occur at Poisson times with rate N, the master equation in weak 
form is, for all test function <f N : Z i —> f N { Z): 




i,j =1 


f N ( Z) f {N) {Z) dZ 



N 


■ ■ Z N ) $(v'i - Vj) dv'i 


N 


■ ■ ■ i Xj , Vj ,. 


(j) N { Z) f {N \Z) dZ. 


( 1 ) 


By exchanging the notations iq and v[ we obtain the following master equation in 
the strong form: 


N 


d,f m ( Z) = Y /(")(Z) (v, • V,,) + JVL/W(Z) , 


i=1 


where the operator L is dehned by 



i,j =1 

*7 


Lemma 4 (Invariance under permutation). Define for all a G &n, 




Then we have: 


L(af^)=a(Lf^) 
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As a consequence, if f^ N \t) is permutation invariant at timet = 0, i.e. crf( N \t) | t=0 = 
f {N) (t) |t =0 for all a G & N , then it is permutation invariant for all times. 


Proof. To emphasise the dependence in Z, we can rewrite the operator L as: 


A 


Lf m ( z) = -jj Y, <( S ( Z ))W( Z ) - v$(Z))a/ w (z) - f {N \ Z), 


N 

*0 = 1 
* 7 H 

with x(Z) = x, Vi(Z) = Vi and 

First note that, setting a Z = (^cr(i), • • •, Z^ N )), we have 

Vi(Z) = V a (i)(aZ) , and P t (af^)(Z) = P a(l) f^(aZ) . 

Therefore by applying the a -1 permutation to the double sum and using the permu¬ 
tation invariance of the rank, see Lemma [2] (ii), we obtain 

1 N 

L°f m ( z> = - <(x< z))i(vi(z) - u(z» fiG/ (,v, )(z) 


*,j=i 


A 


- (<7/ (W) )(Z) 


- V K*V(i, j)](x(Z)) i(V5(Z) - U(Z)) P a(i) f N \aZ) 


* 0=1 


- (<*/ (,v) )(z) 


A 


IV 


E *r A '[r w (<7- 1 (<'),ff- , 0'))](x(Z))i(^-. (j . ) (Z) 


*W 


A 


- U-.(«(Z))* , i'/ (W) (<*Z) - (<t/(">)(Z) 


- V A-"[r' v (i',/)](x(aZ))i(U(aZ) 


*'o'=l 

*w 


- VjfaZ^Pi,f {N \aZ) - f {N) (aZ) 
= (LfW)(aZ) = aLd t f( N \Z). 


The above property states that cr -1 Lcr = L, for all a G ©tv. Supposing that L is 
a bounded operator, we deduce that <r _1 L fc cr = L k , for all A; G N and consequently 
o r_1 e I 'cr = e L . Now, the solution of the problem d t f ^ = NLf (N l with f^\t=o = 
fl, N) can be written f^ N \t) = e NLt f^ N \ We deduce that of^ N \t) = e NLt af^ N \ 
Therefore, if ct/q^* = f$ N \ then <jf^ N \t) = f^ N \t), for all t > 0. If L is not bounded, 
the same property remains true thanks to an approximation argument. □ 
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3.2. First marginal equation for indistinguishable particles. In the remain¬ 
der of this article, we will suppose that /'is invariant under permutations which 
physically means that the particles are indistinguishable. This allows us to define 
the A;-particle marginal as 

(Al) fS ) (Z 1 ,...,Z k ,t)= I f {N \Z ir -- ,Z N ,t)dZ k+1 --- d Z N . 

and is still invariant under permutations of (Z 1; ..., Z k ). 

Proposition 5 (First marginal equation for indistinguishable particles). Assume 
dSU). For any test functions satisfying 

(2) f N (Z u --- ,Z N ) = f(Z 1 ), 

we have 


/ /«( z ) (<VV,MZ) dz 

i =1 ■' 

+ (JV - 1) J T^ 2 (x)^(x 1 ,v 2 )f^(Z) dZ 
+ (N-1)J ^(x)<KZi)t m ( Z)dZ 
+ (JV - 1)(JV— 2) Jw&WftZMWlz) dZ 

-N J fS\Zi) dZi. 


equation (JT|) gives 


cases 

i = 

1 * 

3, 3 = and i > 2, j 


N 


N N N 

dZ = 

■Y. 

4 

+5>?’+E E A, 


3= 2 


i=2 i=2 j=2J4 




-N j «, 

Af 

-) 

~4 

tx) ^(arx, w#)/ <JT) (Z) dZ, 

A 2 ’ 

FI 


;x)^(Zi)/ w (Z) dz. 

■A-iJ 

Al 


; x )0(z,)/ m (z) dz. 


with 
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To compute the first term Aj \ we perform the change of variables Z' 2 = Zj and 
Z'- = Z 2 , which leads to: 


A j 1} = f ■ ■■,x 2 ,...,x N ) (j)(x 1 , v' 2 ) f {N) (Z 1: Z' j ,...,Z' 2 ,..., Z N ) 

dZ 1 dZ'j ... d Z' 2 ... d Z N 

Using the permutation invariance of the rank (see Lemma [2] (ii)), we have 
7T^(xi,x', . . .,X 2 , . . . ,X N ) = 7T%(X i,x' 2 , ...,Xj,.. .,x N ). 


Therefore, dropping the primes and using the permutation invariance of f^ N \ we 
obtain 

Af = J 7r^(x)0(xi,n 2 )/ (JV) (Z) dZ, 

which does not depend on j. 

Similarly, we have 




• ,x 2 , ■ ■ .,x N ) 0(Zi) f [N \Z 1 , Zi,...,Z 2 ,..., Z N ) 

dZ\ dZi... dZ 2 ... dZjsr 


so that, using the permutation invariance of the rank, see Lemma [2] (ii), and the 
permutation invariance of as previously, we obtain 

Af 1 = yb|;(x).XZ 1 )/ < ' v) (Z)dZ, 


which does not depend on i. 

Also, we have with i > 2, j > 2 and i ^ j: 


Aij = j < (x uXi , x j ...,x 2 ,...,x 3 ,...,x N ) 0(ZO 

f^(Z 1 ,Z i ,Z j ,...,Z 2 ,...,Z 3 ,...Z N ) dZ, dZ l dZ J ...dZ 2 ...dZ 3 ...dZ N . 

Then using the permutation invariance of the rank, see Lemma H (ii), and the per¬ 
mutation invariance of as previously, we obtain 

A,j= |ir 2 w 3 (x)^(Zi)/ (A,, (Z)dZ. 

For the last term of (j3J) we obviously have 

J ^(ZiJ/WfZ) dZ = J ^Z t ) fS\Zr) dZL 

Collecting all these identities, we obtain the identity stated in Proposition [5] □ 
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3.3. Propagation of chaos. Assume now that the propagation of chaos property 
holds true he. 

N 

(A2) fW(Z ir -- ,Z N ,t) =n fS\Ze,t), VZ e R 2nN , Vt G [0, oo) , 

i=i 


and define: 



We remark that ply is a probability density. 
We have the following proposition: 


Proposition 6 (First marginal equation with propagation of chaos). Assume (1A2[) . 
For any test functions satisfying (|2]) , we have 


W d t lfS ) (Z 1 )^(Z 1 )dZ l 
» r 

= E / A (dz + (a n ) + (b n ) + (c x ) + (d n ), 


i=l 


with 


(A n ) = 


S N (K) 


<Kx u v 2 )fS\z x )f%\z 2 ) 


N 


K (r N (l, 2)(x)) ]^[ pll } (xi) dx e dZi dZ 2 , 

e=3 

i r N 

( B N ) = ffNjjq J <K z i) fS\z i)K (^(2, l)(x)) (*<) dZi, 

yy or n 

(C N ) = J 0(Zi) (r"(2,3)(x)) IIpS’M dZ t , 


(£>") = -JV J 0(Z t ) /"(Z0 dZi, 

where S N (K ) is given by 


N -1 


S (K) = 


N — 1 


E k 

k =1 


AT- 1 


We note that S N (K ) is the Riemann sum approximation of [J Jl (r) dr. Since we 
assume [J K(r ) dr = 1, S N (K) converges to 1 as iV goes to oo. 

Proof. This result is a direct consequence of Proposition 0 integrating in v when 
possible. We then use that 

K 


K n = 


(N-1)S n (K) ’ 


to obtain the stated result. 


□ 
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4. Limit equation 

For a density p, define the partial mass of p centred in x and of radius s by: 

M p (x, s)= p(x') dx' . 

J\x— af'|<s 

We now state the main theorem of this article: 

Theorem 7 (Limit equation). Assume (1 A2 ji . If 

lim f$ -A f and lim p$ -A p = f f dv, 

N—>00 N-t-OC J 

then, in the limit N —y oo, for all test functions <fi we have: 

d t J /(Z)0(Z)dZ 

= J (t>(xi,v 2 ) f{Z 2 )p(x 1 ) K (M p (xi, \x 2 - Xi|) dxi dZ 2 

- J f(z,) f(z 0 dz 1 , 

or, in strong form: 

-fijf(x,v) +v ■ V x f(x,v) = p(x) J f(x',v)K(Mp(x,\x'-x\)) dx' - f{x, v). 

This result will be obtained by passing to the limit when N —> oo in (j3J). To pass 
to the limit in the transport term of (j3J) is classical and we refer the reader to classical 
textbooks on the subject. We divide the proof of this theorem in two sections. The 
first section deals with the two first terms (A N ) and (B N ) of (j4j) while the second 
will deal with the last two terms (C N ) and (D N ) of 03]). To do so we will be using 
the Bernstein polynomial approximation of functions which is a follows: 

Proposition 8 (Bernstein polynomial approximation, [18]). Let f be a function 
defined on [0,1]. The n-th Bernstein polynomial associated with f is defined by 

:=£/(£) 

If f E C 2 [ 0,1] then 

B n {f 5 x) = f(x) + - -f"(x) + o . 

4.1. Evaluation of (A N ) and (B N ). 

Proposition 9 (Evaluation of (A^) and (B N )). Under the assumptions of Theo¬ 
rem [?] we have for large N 

S n {K ) x (A N ) = 

J f(xi,v 2 ) p!i ] (^i) fS\ Z 2 ) K ^M p ( i)(xi, |xi - ar 2 |)) d^i dZ 2 + o(l), 
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and 


s N (iq x (b n ) = 

j (j)(Zi) /^ } (Zi) K (^M p wi x 2 , \xi - £ 2 |)) Pn\x 2 ) dZi dx 2 + o(l). 
To prove this result we first prove the following: 


Lemma 10. Under the assumptions of Theorem [?j we have for N large, 


N 


K (r JV (l,2)(x)) JJpSv } (^) dx e = K ( M p w{ x u \ x i - ^D) + o(l), 


1 =3 


and 


r N 

/ K (r N ( 2,1) (x)) JJ £$(£/) dx t 
^ £=2 

= j K (M p w(x 2 , \x 1 - x 2 |)) Pn\ x 2 ) dx 2 + o(l) . 

Proof. We first give a combinatorial interpretation of the rank and then use it to 
interpret the terms of the statement as expectation. 

• Let us fix x\ and x 2 . The rank r JV (l, 2) is equal to the number of points £ 3 ,..., xn 
belonging to the ball B = B(xi,\x 2 — x^\) = {x : \x — X\\ < \x 2 — X\\} plus one unit, 
scaled by the factor N — 1, i.e. 

r N d 9 )/ x \ _ Ml E ''' ,N ^ : x i e B } Tl 
y > N- 1 

Denote Pr be the probability such that R N ( 1,2) = R where R N ( 1,2) = (N — 
1) r 7 V (l, 2). To have R N ( 1, 2) = R , we have to choose R — 1 particles amongst N — 2 
to lie in B. The probability that one of the R — 1 particles belongs to B is equal to 

P ■= M m(xi, \x 2 - £i|) , 

PN 

while the probability that one of the N — 2 — (i? — 1) remaining particles lies in M, n \B 
is 1 — p. Therefore, 

( 5 ) 


Now, £1 and x 2 being fixed, the quantity 


N 


/ K (r JV (l,2)(x)) JJpSi } (££) d£ £ , 

J e=3 

can be interpreted as the expectation of K (r JV (l, 2)) when N — 2 points £ 3 ,..., xn 
are drawn according to independent identically distributed probabilities with density 
pi\x) d£. It will be denoted E{K (r JV (l, 2)(x))}. 
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By (ED , we compute 


E { A' (r w (l, 2)(x))} = y; K (-UL) (r _ i) P*‘‘ t 1 “ 


M 


Y. k 

R =0 


A+ 1 

M + l 


M 

R 


P R (1 - p) 


M-R 


with M = N - 2. Since, for N large, K ((R + 1)/(M + 1)) = K (R/M) + o(l) 
(remarking that R/M < 1), 

M 


E{A(r"(l,2)(x))}=^A(- 


R =0 


R\ M 


R 


P R (i- P r-"+o(i). 


M-fl 


Using Bernstein’s approximation, Proposition El we obtain 

E{K(r Af (l,2)(x))} = A» + o(l). 

Which is the first statement. 

• The identity 

r N 

/ ^ (^(2, l)(x)) JJpSv } (^) dx e = K (,M ii) (xi, \x 2 - Xi|)) + o(l), 


1=3 


is obtained in an analogous way by exchanging the role of 1 and 2. We then have to 
integrate by p^\x 2 ) dx 2 to obtain the stated result. □ 

of Proposition 0 Inserting the expressions of Le mm afTOlin (A^) and (B N ) we readily 
obtain the stated result. □ 

4.2. Evaluation of (C N ) + ( D N ). 

Proposition 11 (Evaluation of ( C N ) + (D N )). Under the assumptions of Theorem [?J 
we have 

(C N ) + (D N ) = - U(Z,) /P(Z0 dZ, 

- J 4>{Zi) p { n{x 2 ) K \xi - x 2 \) dx 2 dZi+o(l). 

Like in the previous section, to evaluate (C N ) we first transform the term in 
parenthesis: 

Lemma 12. Under the assumptions of Theorem [7| we have for large N 


N 


I< (r N ( 2,3)(x)) JJpSJ (^) 

e=4 

K'(p 23 ) 


= K(p 23 ) - 

1 

+ iV 


N 

>23(1 -P23) 


(l XB( X2,\x2—X3\)( X l)') 

K"(p 23 ) + 2(1 — p 23 )K' (p 23 ) 


+ 0| A 
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where P 23 = M (i) (x 2 , |x 2 — X 3 I) only depends on X 2 and X 3 and 
Pn 

, v _ f 1 ifx 16 B(x 2 , \x 2 - X 3 |) 

X 13 ( 0 : 2 , 1*2 o: 3 |)l 1 ) | 0 otherwise. 

Proof. Similarly to the proof of Lemma [TOl we interpret the quantity 

r N 

/ K (fN (2,3) ( x )) dx t , 

^ 1 =4 

as the expectation of K (r Ar ( 2, 3)(x)) when the N — 4 points X 4 ,..., xn are drawn 
according to independent identically distributed probabilities with density (x) da;. 
Two cases have to be distinguished: 

• First case: if x\ e B(x 2 , |x 2 — x 3 |) — Like in the proof of Lemma[TU] we have 

n ,o ox #0' e I 4 • • •, N} | : Xj 6 B{x 2 , \x 2 - x 3 |)} + 2 
r iV — 1 

Hence, setting p 2 3 = M m{x 2 , \x 2 - x 3 |) =: p, 


N -1 


E{A-(r w (2,3)(x))} = ^ft' 

R=2 

M 

= E /f 


R 


R =0 


N — 1 

R + 2 
M + 2 


N - 3 


R — 2 

M 
R 


pR-2 _ p) N -3-(R-2) 
p R (1 _ 


with M = N — 3. By expanding K , we have, uniformly with respect to R G 
{0, • • • , M} 


ih 


i? + 2 
M + 2 



2 f M — R\ . f R\ 
+ M\M + 2) \Mj 



Since (M — R)/(M + 2) = 1 — M/R + o(l) we obtain uniformly with respect to 

Re { 0,• • • ,M} 


K 


R + 2 
M + 2 



So, we obtain 


E { A (r"(2,3)(x» } = E K (§) EL (1 - 

_R=0 ' / \ / 



(")/d- P r + o (1 ,. 
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Using Bernstein’s approximation, Proposition |8l to K and to p ha (1 — p)K'(p) we 
obtain 


E {A' (r' v (2,3)(x))} = Kip) + LL/> A"(p) + ^LjU-K'lp) + o(1) 


2 M 


M 


( 6 ) 


= Kip) + U—^-K\p) + dLP) A -"(p) + o(l). 


N 2N 

Second case: if x\ fz B(x 2 , \x 2 — ^ 3 )) — In this case, 


r*(2,3) = 


#{j e {4. •., N} : Xj e B(x 2 , |x 2 -x 3 |)} + 1 


AT- 1 


Following the same step as before we compute, with p = M (i)(a; 2 , \x 2 — a^l) 

Pn 

E{A(r w (2,3)(x))} = A (^Z j) P*"’ (1 - p) A '- 3 -' R - 1 ' , 


which we rewrite 


M 


E{(r w (2,3)(x))} = ^A 


R =0 


R+l 


M + 2 


M 

P 


(1 - p) 


M-R 


with M = N — 3. By expanding K, we have 

R+l 

- I = K I 

M 


K 

So, we have 


M + 2 




M 


E{(r w (2,3)(x))}=^A 


A=0 


M 


R\ (M 


M 


M 


M 


M \R 


M 


M 


R =0 


E b- 


2 R 


p R (l-p) M ~ R 

k ,(r\(m 


M \R 


pR { 1 _p ) M-u + o il) 


M-R 


Using Bernstein’s approximation (see Proposition EJ), we obtain 


E{(r K (2,3)(x))} = Kip) + LL/>A"(p) + U^A'(p) + o(l) 


2M 


M 


(7) 


= A-(p) + hh^A'(p) - UM + LEE A>) + „ (1) . 


N 


N 


2N 


• We obtain the result stated in Lemma [12] by noticing that Expression flTj) is equal 
to the sum of Expression ([ 6 ]) and an extra term —K'(p)/N. □ 

Lemma 13 (Evaluation of S N (Ii ) x ( C n )/(N — 2)). Under the assumptions of 
Theorem [?], we have 

m x(c «)W 1 + l + A ' (1) - A W 

N - 2 K ’ V N 2N 


4>{Zi) A 1 } (^i) dZ 1 


Jy j <l>{Zi) Jn\Zi) Pn\x 2 ) K (^M p w ( x 2, \ x 2 ~ ^3^ dx 2 dZ 1 + o(l). 
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Proof. Using Lemma [T2] and separating the cases x\ G B{x 2 , |^2 — ^ 3 1) and x\ fL 
B(x 2 , \x 2 — m 3 1), we can write 

x (c N ) = (1) + (2) + o (1), 

where, writing p for p 2 3 , i.e. p = M (i)(x 2 , \x 2 — x 3 |), we have: 


( 1 ) = Tj I <K^i) fN\ z i) K '(p) P { nM Pn\ x s) da; 2 dx 3 dZi 

>Ia:iGB(a:2,p2— rc31) 



and 


(2) = J Mzj (a-(p) + + hi_d A -» (p) ^ 

Pn\ x 2) Pn\x 3 ) dx2 dx 3 dZi + o 



• For the term (1), we first notice that X\ G B{x 2 , \x 2 — x 3 |) is equivalent to saying 

that X3 4 B{x 2, laq — x 2 \) so that, for p = M (p (x 2 , \x 2 — x 3 |), 

Pn 


« = V 


X3^B(X2,\xi-X2\) 


f^\Zi)K\p) Px\x 2 ) p$(x. 3 ) dx 2 dx 3 dZ 1 




By the change of variable stated in Lemma fl4l and applied to H = K', p = , 

x = x 2 , and r = \x\ — x 2 \ we have 


' X3^B(X2,\X1—X2\) 


K' i)(x 2 , \ x 2 - Z 3 I)) Pn\ x 3 ) d^3 


= K( 1) - K (^M p w (x 2 , \ x i ~ x 2 \} 


Inserting this in (1) we obtain 


(8) N x (1) = K( 1) J <f>(Zi) f%\z x ) dZ 1 j p$(x 2 ) dx 2 

- J 4 >{Zi) fx\Zi) Px\x 2 ) K ^M p (i)(x 2 , \ x i - x 2 \j dx 2 dZi + o(l) 

= K( 1) j tl>(Z 1 )fS ) (Z 1 )dZ 1 

- j 4>(Zi) f^\Zi) Pn\x 2 ) K (^M p w (x 2 , \x\ — £ 2 ^ dx 2 dZ\ + o(l). 


hi), 
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• Using again the change of variable result of Lemma fl4l together with integration 
by parts, we compute: 


K(P) + i-^iC(p) + dL_P> K"(p))pW(x 3 ) dz 3 


2 N 


= / (K(p) + 


1 - 2/1 TS<(~\ , P( l ~P) 


N 


-K'(p ) + 


2N 


K"{p)) dp 

1 K{ 0) + A'(1) 


= 1 + 


N 


2N 


where p = M m(x 2 , \x 2 — xd). From this and Lemma H2l we deduce: 

p N 


1 _ K (0) + K(l) 
N ~ 


(9) (2) = ( 1 + - - ) / <\>(Zf) f^\z x ) dZ x + o 


2N 


• Combining the two terms (JH]) and ([9]), we obtain the result stated in Lemma fl3l □ 
We are now ready to prove Proposition |TT] 

of Proposition [HI The proof is divided in two main steps. 

• We first have 


N—l 


s n (k) = — y k 

y > At — 1 


k 


N - 1 


iV — 1 

fc=i 

K(l)-K(d) 1 

2(N-l) N- 1 


AT(0) + K(l) 


N—2 

Y. k 


k =1 


k 


N — l 


In the second term of this expression, we recognize the approximation of K ( s ) ds 

by the trapezoidal rule. As the trapezoidal rule is second order, it leads to 


2N I N 


As a consequence 


N-2 
S N (K ) 


= N 


1-2/N 


= N - 


1 + (K(l) - K(0))/2N + o(l/N) 
K(l) - K (0) 


( 10 ) 


2 


2 + o(l) . 
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• Now collecting the estimate of Corollary [13] and (HU we obtain 
N — 9 

(CO =5^ Id) + (2)1 

= (jv - /t(1> ~ Jt(0> -2) / <A(Zi) f!v\Zi) d Z 1 

+ (i + /t(1) ~ /f(0> ) J4>(Zi) ftt’izo d z 1 

- f <j>(Z 1 )f$(Zi)p§\x 2 )K(M p cL)(x 2 ,\x 1 -X 2 \') dx 2 dZ 1 

+ o(l) . 

And as 

(. D n ) = -N f ftZj)dZ l7 

we obtain the statement of Proposition dH □ 

4.3. Proof of Theorem 0. We have to pass to the limit in (j3J) . By Propositions [9] 
and dn we have 

r N . 

St / /^(Zi)^(Zi) d Zt-J2 fS\ Z) («<' V,,)<#.(Z) dZ 
= ( +N(/v) _ X )/ K l x i _a: 2 |)) pS } (x 2 ) dZi dx 2 

+ J 0(xi,T 2 )p!i ) (^i) A 1} (^ 2 ) AT (a^(i)(®i, |a?i—a? 2 |)) dD dZ 2 

- J tl,(Z 1 )fS\z 1 )dZ 1 + o(l). 

As N goes to oo, the second line goes to 0 since S N (K) is the Riemann sum approx¬ 
imation of Jq K(r) dr = 1. The convergence in the other terms is formally obvious 
and leads to the stated result. 


5. Discussion 

5.1. Large-time behaviour. Consider a function homogeneous in space (t,x,v) i— >• 
G(t,v). Since f K = 1, by Lemma dH we get 

f)C f 

—(v) = -v ■ V x G(v) + G(v) / K ( M p (x , \x' - x|)) dx' - G(v) = 0 . 

Hence any function homogeneous in space (t,x,v) H > G(t,v ) is a stationary solution. 
Moreover, on a periodic spatial domain, we can expect that any solution converges 
at large-times toward a function of this type. The proof of such a claim is left to 
future work. 
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5.2. Discrete versus continuous approach. We can wonder if the large-time and 
large number of particles limits permute. It does not seem the case. Indeed, the 
number of distinct velocities decreases when there is a finite number of particles 
while, as discussed in the previous section, the distribution of velocities remains 
constant in time in the case of a continuum of particles. 

In the case of a finite number of particles the consensus in the direction the particles 
adopt is longer and longer to obtain, see Figures [T] [2]and|3j 




FIGURE 1. Trajectories of the particles on the left, variance and number of 
different speed as functions of time in the case of 10 particles taken randomly 
in [—10,10] for the position and for the speed. 




FIGURE 2. Trajectories of the particles on the left, variance and number of 
different speed as functions of time in the case of 20 particles taken randomly 
in [—10,10] for the position and for the speed. 


6. Conclusion 

In this paper, we have investigated a system of particles interacting through leader 
following interactions where the choice of the leader is determined by a topological 
rule. Under a propagation of chaos assumption, we have shown that the large system 
size limit is described by a spatially nonlocal kinetic model of Boltzmann type. This 
result heavily relies on approximation properties of Bernstein polynomials. Obvi¬ 
ously, the very simple leader following model considered in this paper offers many 
directions of complexihcation leading to biologically or socially more realistic rules. 
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FIGURE 3. Trajectories of the particles on the left, variance and number of 
different speed as functions of time in the case of 70 particles taken randomly 
in [—10,10] for the position and for the speed. 


An example could be the introduction of some noise, e.g. the velocity after the in¬ 
teraction would be randomly selected according to a probability law centred around 
the leader velocity. One could also think of the two particles joining their average 
velocity up to some noise, in the spirit of [2]. Finally, binary interactions with the 
closest neighbour could also be investigated. 


Appendix A. Fundamental lemma 
Lemma 14. For any H, 

r i‘Mp(x,r) 

I(x;r) := / H(M p (x,\x'— x\)) p(x') dx'= / H{p) dp . 

J B r (x) Jo 

Proof. First note that since 

M p (x,s) = 

we have 


s<s Ju jsS ™- 1 


p(x + su) s 11 1 ds da;, 


d f 

——M p (x,s)= / p(x + so;) s n_1 da;. 

dS gn-i 


Using the polar change of variables, 


\x — x =: s 


x — X 

] / r ■ ^ > 

\x' — x\ 


we have 


I(x;r) = / / p(x + su) ff(M p (x, s)) s n 1 ds do; 

Js<r J c<jE § n—1 

= [ —j-M p (x, s) H(M p (x, s)) ds . 


f s<r 


ds 


Setting p = M p (x , s), so that dp = - 4-M p (x , s) ds, we obtain the stated result. □ 
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